Critical Behavior of Light 
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Light is shown to exhibit critical and tricritical behavior in passive mode-locked lasers with ex- 
ternally injected pulses. It is a first and unique example of critical phenomena in a one-dimensional 
many body light-mode system. The phase diagrams consist of regimes with continuous wave, driven 
para-pulses, spontaneous pulses via mode condensation, and heterogeneous pulses, separated by 
phase transition lines which terminate with critical or tricritical points. Enhanced nongaussian fluc- 
tuations and collective dynamics are observed at the critical and tricritical points, showing a mode 
system analog of the critical opalescence phenomenon. The critical exponents are calculated and 
shown to comply with the mean field theory, which is rigorous in the light system. 

PACS numbers: 42.55.Ah, 42.65.-k, 05.70.Fh 



Statistical light-mode dynamics (SLD) 0, 0] that is 
based on statistical mechanics, provides a powerful ap- 
proach for the study of complex nonlinear light systems, 
while serving as a new statistical mechanics paradigm. It 
was developed to treat and solve long standing questions 
in laser physics. A prime example is the laser pulsa- 
tion threshold, a phenomenon in passive mode locking 
that attracted much attention Passive mode- locking 
occurs when a saturable absorber element is placed in 
the cavity, driving the laser to generate extremely short 
light pulses. The transition mechanism from a continu- 
ous wave (cw) to pulsation of such lasers was shown via 
SLD 0, Q to be an intrinsic property of the many inter- 
acting mode system, ruled by the balance between the 
nonlinear interaction induced by the saturable absorber 
and the randomizing effect of noise. The theory, verified 
by experimental study, demonstrated the existence of 
first order phase transitions between disordered (cw) and 
ordered (locked) mode phases, as the noise ("tempera- 
ture") or the laser power were varied. Thus noise alone 
stabilizes the cw state, showing a noise induced phase 
transition. p(. 

In the present Letter we report on findings of critical 
and tricritical phenomena in the many light-mode sys- 
tem. It is a first example of such behavior with light, 
that also provides a physical realization of a strict one- 
dimensional many body system. For criticality to appear 
in the laser mode system we add to it an external driv- 
ing field which is a counterpart of the external magnetic 
field in magnets and the pressure in gas-liquid-solid sys- 
tems [6|. It is achieved by injecting the laser with pulses 
from an external source, which in the simple case matches 
the repetition rate of the laser. When the injection is 
weak, the ordering phase transition persists, shifted to 
higher "temperature", with a "para-pulses" phase (pulses 
induced and driven by the external injection). However, 
the phase transition line terminates in a critical point, 
where the distinction between para-pulses and sponta- 
neous pulses disappears, similarly to the vapor-liquid 
critical point. Thus, by increasing the external injection, 



it is possible to obtain mode locking smoothly from cw. 
Near the critical point the system exhibits the familiar 
critical phenomena, including divergence of response co- 
efficients, characterized by universal critical exponents, 
and nongaussian critical fluctuations, enhanced by a fac- 
tor of N 1 / 4 , where N is the number of active modes, 
compared to normal fluctuations. The latter could mean 
a two orders of magnitude fluctuations enhancement in 
practical systems. It is a light mode system analog of the 
critical opalescence phenomenon. 

When the injection repetition rate is higher than that 
of the laser, the high- and low-temperature phases are 
characterized by equal and unequal pulse powers, respec- 
tively, which cannot be smoothly connected. However, 
beyond a threshold injection level the transition becomes 
continuous rather than first order. The two phase tran- 
sition lines meet at a tricritical point, often found in 
systems with interaction competing with external driv- 
ing jjj, around which tricritical behavior is observed, with 
its distinct set of universal exponenets, and the tricritical 
fluctuations are enhanced by a factor of N 1 ^ 3 . 

The inclusion of external injection has other important 
consequences. For example, as in the vapor-liquid case, 
field induced "condensation", that is in our case mode 
ordering and pulsation, can be sustained to some extent 
when the field is removed, only below the critical point. 
Another practical aspect of external injection is the pos- 
sibility to overcome basic deficiencies in passively mode 
locked lasers, their low repetition rate and timing jitter 
. Our theory provides a basic and quantitative under- 
standing of the way that an external field locks the laser 
pulsation, and the threshold power needed (or maximum 
noise allowed) for various operation scenarios to happen. 

We perform our theoretical study in the framework of 
the coarse-grained model of SLD, representing the cavity 
electric field envelope ip by a single variable in an inter- 
val whose length is of the order of the pulse width. The 
derivation of the model from the passive mode locking 
master equation has been discussed before 0, 0, [1 , 
and it shows good quantitative agreement with theoreti- 
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cal studies of the master equation as well as with exper- 
iments. 

The dynamics of the the field variable in interval num- 
ber n is expressible as 



dH 

dip™ 



+ gipm + r) m (t) ■ 



(1) 



The "Hamiltonian" H with injection strength h takes the 
form 



JY 



H =-iVRe^(|l^m| 4 + 2/ lm C) 



(2) 



where N 3> 1 is the number of active modes, 7 is the 
coefficient of saturable absorption, and h n is the external 
injection at site n. h n will be assumed to take nonzero 
values only at a small number n of intervals, g is the over- 
all net gain, which can be assumed without loss of gener- 
ality 0] to set the intracavity power ||V>|| 2 = J2 n IVvJ 2 to 
a fixed value P, in which case g becomes a Lagrange mul- 
tiplier for the constraint. The random term 77, represent- 
ing noise from spontaneous emission and other sources, 
is modelled by a (complex) Gaussian white noise with 
covariance «(i)7?m(*')) = 2T5 nm S(t - t'). 

The invariant measure of Eq. (01 is a Gibbs equilib- 
rium distribution [HHll| 



p[i) ] = le-W8{U\\ 2 -P), 
with the partition function 

Z - / [d<ip][dip*]e- H M/ T 5(\\i>\\ 2 - P) 



(3) 



(4) 



As in previously studied cases of SLD [1|, |2j, |4| , when 
N ^> 1 the invariant measure is concentrated on config- 
urations where all but a finite number of the tp variables 
are 0{N~ 1 / 2 ). Here the intervals where -0 = 0(1) are 
precisely the n intervals which are subject to external in- 
jection. The free energy F — —TlogZ is then given by 
F = ATmin^,...^ /„(V>i, ■ ■ - ,ipn), where 



fn = - ^(^|V m | 4 +2Re/ lm V™)+Tlog(P-^ |^ 



rn — 1 



m — 1 



(5) 

We establish Eq. (J3J) using the results of |2( for the par- 
tition function Z for the h = case, 



Z 



d 2^ e [N 7 W 4 +T\og(P-W 2 )]/l 



(6) 



asymptotically for large N. We can proceed immediately 
to perform the integration in Eq. (0} over the N — n 
intervals where h = 

Z~ f d 2 ^]Jd 2 ^ m (7) 

xe -T(.H(4, 1 ,...,ip n ) + N 1 \i,\ 4 +NTlog(P-Y;^ =1 \4> m \ 2 -M 2 ) 



The exponent in the integrand in Eq. J7J) is proportional 
to the large parameter N, and therefore the integration is 
concentrated near the global minimum of the integrand. 
It is straightforward to verify that the minimum is al- 
ways obtained when ip = 0, which then implies that 
F = Nminf n . 

Moments of the pulse strength can be obtained in the 
standard manner by taking derivatives of the free energy. 
Alternatively we may obtain directly the probability dis- 
tribution function of the pulse amplitudes by integrating 
out the N — n unforced variables in Eq. © getting 



P(^,...,^ n )~e^-^--*») 



(8) 



The width of the distribution tends to zero in the ther- 
modynamic limit N — » 00, which shows that the pulse 
amplitudes are thermodynamic observables. 

We henceforth specialize to the case that all nonzero 
injection values arc of the same magnitude h, that is 
appropriate for injection of pulses from a source with a 
repetition rate n times faster than that of the laser. For 
the minimization problem there is no loss of generality in 
assuming that the injection values are all real, since the 
minimizing ip values have the same phase as the corre- 
sponding h values. By a simple rescaling of the variables 
the free energy /„ can be reduced to 



/ 



n 



t)X r , 



Tlog 1 



22 x m 



(9) 



where real and positive. It follows from Eq. © 

that thermodynamics depends on only two dimension- 
less parameters, the reduced temperature T = -^p? (the 
inverse of of the interaction strength |l|) and reduced 
driving [) = p^ 2 . x m , the minimizers of /, are related 



to the expectation values (in the invariant measure) of 
the pulse powers by x 2 ^ = (\ip m \ 2 ) /P- 

The study of thermodynamics and critical behavior is 
now reduced to the analysis of the function / and its 
minima. We note that for any values of the parameters, 
minima can occur only at configurations where at least 
n — 1 of the equal, with value y, and the other 

minimizer, which we denote by x is greater than or equal 
to y. Accordingly the free energy is obtained from 



f{x,y) 



1 



(rc-l)y 4 ] -t)[x + (n-l)y] 

-Tlog[l-x 2 ~(n-l)y 2 llQ) 



It is instructive to consider the above results in the k 
or mode space a/c, the discrete Fourier transform of ip n . 
The Hamiltonian reads 



H 



a i a t a l a l 

j — k-\-l — m— 



2ReJ2 h k a* k , (11) 



where hk takes nonzero value for k being integer multiples 
of n. Therefore, when n > 1, the set of modes consists 
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of two types, with and without the presence of the ex- 
ternal field. The hamiltonian in Eq. is analogous 
to the one of an antiferromagnet placed in an external 
homogeneous magnetic field |8| , but with the roles of the 
interaction term and the driving term reversed. Namely, 
the driving acts on a subset of the modes, but the in- 
teraction tends to align all modes in the same amplitude 
and phase. The consequences of this competition are de- 
rived below. The free energy can be expressed in the 
mode representation using the mode amplitude expecta- 
tion values (a/) and (a u ) of those with (forced) and with- 
out electric field (unforced), respectively. The relations 
are x = £((a/) + (ra - 1) (a u )) and y = £((a/) - 1 (a u )) 
and the free energy is 

f = _ i [(a/)4 + {n ~ 1)(n2 " 3n + 3) K)4 

+6(n - 1) (a f f (a u f + 4(n - l)(n - 2) (a/) (a„) 3 ] 
-[) (a f ) - T log (1 - 1 (a f f - ^ (a u f) (12) 

For the analysis below we use the real space formula- 
tion. We consider first the case of n — 1, where there is 
a single pulse (the external field is applied on all modes), 
and / (Eq. 10) depends on the single variable x. For 
zero and small values of f), / has a minimum X\ near zero 
and, for small enough T another minimum X2 > x\ be- 
low 1. For such () there is a threshold temperature Ti(t)) 
where f(x%) = f(x2). As T is decreased through this 
line x jumps from x% to xi in a first order phase transi- 
tion. When t) > the jump is between two pulsed states, 
but the high-temperature phase pulses are driven pulses 
whose power decreases smoothly to zero when \) — > 0. 
We term this phase "para-pulse" because of its resem- 
blance with the paramagnetic phase of magnets above 
the Curie point. For sufficiently strong t), on the other 
hand, / always has a single minimum x between zero and 
one, which decreases smoothly from one to zero as T is 
increased. Thus, the coexistence line 71(f)) terminates at 
a critical point (t) c ,T c ). The phase diagram is shown in 
the upper part of Fig. ^ 

The lower part of Fig. ^ shows function of T 

for several \) values, x undergoes a jump for f) < f) c , and 
displays an infinite slope at the critical point — a mani- 
festation of the critical divergence of the susceptibility. 
The critical point itself is characterized by the vanishing 
of the first three derivatives of /, which gives three poly- 
nomial equations for the three unknowns x c , f) c and T c . 
The equations can be solved explicitly by radicals giving 
x c re 0.53, f) c « 0.2035, T c w 0.34. 

One may define the usual critical exponents [(| j3, 7, 
5 by (x - x c )| cooxis tcnco ~ (T c - Tf, \ = cfx/dt) ~ \T - 
Tc]" 7 , an d (x -x c )\t=t c ~ (f) - fk) 1 / 5 - The exponents, 
as well as the nonuniversal amplitudes can be calculated 
by the standard procedure of expanding / up to third 
order near the critical point [6| , which yields the classical 
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FIG. 1: Upper figure: The phase diagram n=l (homoge- 
neous injection to all modes that matches the cavity repe- 
tition rates). The lines gives the first order phase transition 
curve, that terminate at the critical points. The optical phase 
4> follows the phase of f). Lower figure: Pulse power vs. T for 
different values of external injection f). 

mean field exponents /3 = 1/2, 7 = 1, S = 3, expectedly, 
since mean field theory applies to our system. 

The fluctuation-dissipation relations naturally hold in 
the laser mode system; it follows from Eq. 10} that 

that is, the critical exponent 7 also describes the diver- 
gence of pulse power fluctuations near the critical point. 
To study the fluctuations at the critical point we turn 
back to Eq. JHJ, which in the present context is a proba- 
bility distribution for the single pulse amplitude. Letting 
x = P~ 1 l' 2 \ij)\, the criticality condition implies that for 
T = T c , h = h c 

P(x) ~ e^K^c^+otz-zc) 5 ) ^ (14) 

where a is O(l). The fluctuations distribution is non- 
gaussian, and the scale of the critical fluctuations is 
0(iV -1 / 4 ), stronger by a factor of TV 1 / 4 than normal 
fluctuations. The critical fluctuations are much larger 
than the typical amplitude of the continuum background, 
0(iV -1 / 2 ). It follows that the fluctuations of the contin- 
uum background are correlated, being the SLD analog of 
the critical opalescence phenomenon. 

The thermodynamics with n > 1 is qualitatively differ- 
ent. The external injection encourages the formation of n 
equal pulses, clashing with the tendency of the saturable 
absorber to form a single strong pulse. As a result, the 



phase diagram consists of unequal pulse phase for weak 
noise and weak injection, and an equal pulse phase for 
strong noise or strong injection. As the two phases are 
characterized by different symmetries, there can be no 
smooth transition between them, and they are separated 
by a phase transition line. However, the phase transition 
may be continuous or first order, depending on whether 
the transition order parameter q = (x — y)J y/x 2 + y 2 is 
continuous or jumps to a nonzero value at the transition. 

We consider in detail the case n — 2 where both be- 
haviors occur. It is more convenient to express / of Eq. 
(I10|l in terms of q and p = x 2 + y 2 



f(p,q) = -^(l+2g 2 -g 4 )-h^2p(l - ^)-Tlog(l-p) 

(15) 

to be minimized over p and q. The values of the minimiz- 
ers p and q, giving their "thermal" averages vs. T, are 
given in Fig. 2. / is manifestly symmetric in q, since the 
pulse amplitudes x and y play symmetric roles, where- 
from it follows that / is always stationary with respect 
to q when q = 0. The condition d q f = has another so- 
lution p 3 = \{q 2 — 1) 2 (2 — q 2 ), and the global minimum 
of / is reached in one of these configurations, the first 
corresponding to equal and the second to unequal pulse 
mode locking. 

Straightforward analysis shows that for large T the 
function / has a single minimum, which occurs at q = 0. 
For large f) this situation persists as T is lowered until at 
T = % = |fj 2 / 3 (2 — f) 2 / 3 ) the minimum becomes a saddle 
and two minima with nonzero q form, i.e., q undergoes a 
continuous phase transition, see Fig. |3| For small t), on 
the other hand, nonzero q minima appear for T > 7j,(fj), 
and at Ti(\j) exchange stability with the q = minimum 
in a first order phase transition, also shown in Fig. |3| 

The intersection of the line of first order phase tran- 
sition 7i(fy) and the line of continuous phase transition 
7fc(fj) can be shown to occur at (7^,f) t ) = (^jg-,^)- 
(T t ,\)t) is a tricritical point 0, with symmetric tricrit- 
ical phenomena, and the phase diagram Fig. [21 is quite 
similar to that of metamagnets ||, where tricritical be- 
havior is known to occur. In particular we may define 
tricritical exponents such as fit and fct associated with 
nonsymmetric (e.g. q) and symmetric (e.g. p) fields re- 
spectively by q ~ {% - T) ft and p - p t ~ (% - T) tht 
near the tricritical point, see Fig. [5] As before, the ex- 
ponents take the classical values f3 t = 1/4, /3 2 t = 1/2 
pj. Near the continuous phase transition line there are 
ordinary critical phenomena, for example q ~ (% — T) 13 
and p-p b ~ {%- T) fe , where /3 = 1/2 and fc = 1. At 
the tricritical point fluctuations are enhanced by a factor 
jV 1 / 3 compared to normal fluctuations. 

When the ratio of the repetition rates n is three or 
larger one can show that the transition between the equal 
and unequal pulse phases is always first order; a typi- 
cal phase diagram is shown in the right panel of Fig. [3] 
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FIG. 2: (x — y) (left) and x 2 + y 2 (right) vs. the normalized 
temperature for different values of f) for injection repetition 
rate twice the laser repetition rate. The nature of the phase 
transition changes at t) = t)t 



Critical and multicritical phenomena could be observed 
in these cases under external injection of unequal pulses. 
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FIG. 3: The t) — T phase diagrams for repetition rate ratios 
n = 2 (left) and n — 20 (right). The bold line is a first order 
phase transition curve. The dotted and dashed lines form 
together the bifurcation curve of which the first is a 

line of continuous phase transition. 

The continuation of the first order phase transition 
line (in all cases) under external injection to power lev- 
els below the one required for spontaneous passive mode 
locking (self-starting threshold) has an interesting prac- 
tical implication. A laser operating below the threshold 
pumping level may be mode locked by external injection. 
The injection may be removed and, provided pumping is 
strong enough, the laser would remain metastably mode 
locked for an exponentially long lifetime Q. 
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